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Abstract 

Recently, there has been some interest for building supersymmetric models of dou- 
ble inflation. These models, realistic from a particle physics point of view, predict 
a broken-scale-invariant power spectrum of primordial cosmological perturbations, 
that may explain eventual non-trivial features in the present matter power spec- 
trum. In previous works, the primordial spectrum was calculated using analytic 
slow-roll approximations. However, these models involve a fast second-order phase 
transition during inflation, with a stage of spinodal instability, and an interruption 
of slow-roll. For our previous model of double D-term inflation, we simulate nu- 
merically the evolution of quantum fluctuations, taking into account the spinodal 
modes, and we show that the semi-classical approximation can be employed even 
during the transition, due to the presence of a second inflaton field. The primordial 
power spectrum possesses a rich structure, and possibly, a non-Gaussian spike on 
observable scales. 

PACS: 98.80.Cq; 04.62.-hv; 05.70.Fh 
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1 Introduction 



The standard inflationary prediction [1] concerning the primordial power spec- 
trum of scalar metric perturbations (denoted $ in the longitudinal gauge) is 
a simple power-law : < |$p> oc k^~^. However, for more than ten years, 
there has been some interest for models now called broken-scale-invariant 
(BSl), predicting deviations from a power-law. These models generally in- 
volve, in addition to the usual inflaton field, other (effective) fields, driving 
successive stages of inflation, or just triggering phase transitions [2-9]. From 
the beginning, BSI models have been motivated by observations. After decisive 
experiments (like the APM redshift survey and COBE), it appeared that the 
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standard cold dark matter model (sCDM), even with a power-law primordial 
spectrum, was not in agreement with observations. Adding the BSI primordial 
spectrum of "double inflation" could improve the situation [10,11], but later 
it was shown that this model cannot account for small-scale cosmic microwave 
background (CMB) anisotropy measurements [12]. Nevertheless, interest for 
BSI models revived recently, when some authors found possible evidence for 
a feature in the present matter power spectrum, around 120 h~^Mpc [13-16] : 
this is still a controversial and open question, which future redshift surveys 
should answer precisely; anyway, in [17], it was shown that this feature could 
be associated with the BSI step-like spectrum of an inflationary model pro- 
posed by Starobinsky [18] (combined with a cosmological constant), while in 
[15], similar conclusions were reached with a spike in the primordial spectrum. 
Finally, there has been a recent attempt [19] to connect BSI models with a 
possible deviation from gaussianity in COBE data [20]. 

The step-like spectrum of [18], based on a non-analyticity in the inflaton 
potential, intends to be an effective description of some underlying, more 
complicated and more realistic model. As noticed in [21], this underlying 
model should certainly involve more than one scalar fleld, and some rapid 
phenomenon, such as a phase transition, occurring approximately 60 or 50 
e-folds before the end of inflation. Such models were recently proposed [22-24] 
in the context of global supersymmetry or supergravity.Q The main advantage 
of supersymmetric inflationary models is that they do not require very small 
parameters, due to the existence of flat directions in the tree-level potential, 
and of small one-loop quantum corrections in the effective potential. 

In these works [22-24], the power spectrum was calculated using analytic 
formulas, for modes exiting the Hubble radius during one of the two slow- 
roll stages. Indeed, when inflation is driven by a single slow-rolling scalar 
fleld, a simple expression relates the spectrum of scalar metric perturbations, 
< |$P >, to the fleld potential (and its derivatives) at the time of Hubble 
crossing for each mode. As we will see later, in some situations, the analytic 
results can be extended to multiple slow-roll inflation. However, in more com- 
plicated situations, it is necessary to calculate the primordial spectrum from 
flrst principles, i.e., to follow numerically the evolution of the mode functions 
of fleld and metric perturbations : this amounts in integrating a system of 
second-order differential equations. At the end of inflation, the mode func- 
tions of $ give the primordial spectrum. Here, we want to perform such a 
numerical simulation, in order to check the analytic results, and to obtain the 
shape of the primordial spectrum on intermediate scales, for modes exiting 
the Hubble radius between the slow-roll stages. For these scales, and also on 



Also, while this work was being completed, there has been an interesting proposal 
for generating features in the primordial power spectrum through resonant fermion 
production during inflation [25]. 
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larger scales, the calculation appears to be very interesting, and reveals some 
unexpected features. 

The supersymmetric double inflationary models of refs. [22-24] are, in the 
vocabulary of inflation, double hybrid models,[^ involving two inflaton flelds, 
and two "trigger" flelds, driving second-order phase transitions : one between 
the two inflationary stages, and one at the end of the last stage. Since we 
will focus on scales crossing the Hubble length around the intermediate phase 
transition, the second trigger fleld, the one ending inflation, will not be im- 
portant.!^ So, our study involves three flelds : the inflaton fleld A, in slow-roll 
during both stages of inflation; the inflaton B, in slow-roll during the flrst 
stage only, and then performing damped oscillations; flnally, the trigger fleld 
C, stable in a local minimum during the flrst stage, and then performing oscil- 
lations around the true minimum. In fact, C is a complex Higgs fleld, charged 
under a U{1) gauge symmetry, and the phase transition describes the sponta- 
neous symmetry breaking of this U{1). The inflatons A and B are also complex 
flelds, but their phases are not affected by the evolution, and are flxed from 
the beginning; so, we can treat them as real flelds. 

This fast second-order phase transition, triggered by the non-slow-rolling fleld 
C, and taking place approximately 50 e-folds before the end of inflation, while 
the fleld A is still in slow-roll, is a new situation for the calculation of the 
primordial power spectrum. 

Indeed, in the context of inflation, second-order phase transitions have been 
flrst considered by Kofmann and Linde [4] and Kofmann and Pogosyan [6], 
but for models with two flelds instead of three : an inflaton and a trigger fleld. 
Then, during a short stage, corresponding to the beginning of the phase transi- 
tion, the mass of the trigger fleld becomes negative, and adiabatic/isocurvature 
perturbations are exponentially amplifled. This mechanism is often called spin- 
odal instability. It results in the formation of a large and narrow spike in 
the primordial spectrum, together with the appearance of topological defects, 
which are diluted later if the second stage of inflation is sufficiently long [4]. 
These models were recently studied in more details, in the context of "super- 
natural inflation" [29] , and "natural hybrid inflation" [30] , and connected with 
primordial black hole production. However, with two flelds only, there is not 
as much freedom as in our model (in particular, because the trigger fleld must 

^ Double hybrid or double supersymmetric inflation had been introduced earlier 
[26], not in the context of BSI primordial spectra, but as a solution to the problem 
of initial conditions for hybrid inflation. 

^ However, the second phase transition will produce cosmic strings [27], which can 
be as important for the formation of structure as the primordial power spectrum 
itself [27,28]. In current studies, these mechanisms are completely separated, and 
here we consider the calculation of the primordial spectrum only. We are grateful 
to Rachel Jeannerot for pointing out this fact. 



3 



be light enough for driving a second inflationary stage). Then, the spike in 
the primordial spectrum produces enormous density fluctuations {5p/ p ~ 1), 
allowed only on very small scales; so, the second stage of inflation cannot last 
more than 30 e-folds. In our model, since we have a third fleld driving the sec- 
ond stage of inflation, we will observe for natural parameter values a moderate 
spike, which exceeds the approximately flat power spectrum only by a factor 
of order one. This is typically the kind of feature that can improve current 
flts to the observed matter power spectrum [15], and in contrast with previous 
models, it is worth calculating precisely the spike shape and amplitude. 

Second-order phase transitions during inflation have also been recently studied 
in details by Boyanovsky, de Vega, Holman and collaborators (see for instance 
[31-33]). However, in these works, the Higgs fleld is in slow-roll, and drives 
inflation : it is a very interesting situation, involving strongly non-linear evo- 
lution, and requiring some non-perturbative approach. The stage of spinodal 
instability leads to the emergence of a stochastic homogeneous background 
from the coarse-graining of super-Hubble scales, like in stochastic inflation 
[34]. We will see that in our model, such features are avoided, and the usual 
perturbative semi-classical approach can be employed (again, because we have 
an inflaton A supporting inflation during and after the transition). 

Obviously, the reason for which double hybrid inflationary models did not get 
a lot of attention before, is that without supersymmetric motivations, they 
appear as quite "non-minimal", with many free parameters (some of them 
being tuned to unnaturally small values). On the other hand, if one plays with 
superpotentials and tries to build double-inflationary models, double hybrid 
inflation arises as one of the simplest possibilities. One can imagine many 
variants of double supersymmetric inflation : the potential can be of the type 
of double F-term inflation [22], double D-term inflation [23], or mixed; special 
supergravity effects can arise [24]; different gauge groups can be broken [27]; 
different assumptions can be made about the fleld charges. In the following, we 
will consider only our previous double D-term model, and neglect supergravity 
corrections|3 but most of our results would arise in the other variants, for some 
(natural) parameter values. 

Which kind of novelty do we expect from this study? 

• First, since slow-roll conditions are abruptly violated during the transition 
(so that the Hubble parameter H will depend as much on the kinetic energy 
of the oscillating flelds B, C, as on the potential energy of the inflaton fleld 
A), we should devote a special attention to the evolution of super-Hubble 

^ These corrections are important when the fields take very large values, close to 
the Planck mass. In particular, this may be the case at the beginning of inflation 
[35,36,24]. In our model, with the small gauge coupling constants that we will con- 
sider, supergravity corrections are sub-dominant. 
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wavelengths during the intermediate stage. Previous studies [22-24] assume 
an analytic expression for the amplitude of large wavelength perturbations 
(those which exit the Hubble radius during the first stage), of the type < 
|$P> oc V J2i{Vi/Vi)'^ (see the next section for definitions and references), 
where the sum runs over slow-rolling fields : in our case, A and B. Strictly 
speaking, this formula was originally derived and employed in the context 
of multiple slow-roll inflation [2]. We want to test this expression in our 
case; we will see in fact in section 3 that the violent slow-roll interruption 
induces a significant deviation from the analytic prediction (with our choice 
of parameters, by a factor 4.3 for /c^/^ < |$p>^/^). 
• second, interesting physics could emerge from the description of the Higgs 
field dynamics during the transition. At the very beginning of the phase 
transition, the Higgs field has to roll away from an unstable equilibrium 
point. For the purpose of the simulation, one may think that it is sufficient 
to put by hand an arbitrary small value for the Higgs zero-mode, and then 
let this zero-mode roll down. Actually this is completely misleading. The 
final result would depend strongly on the initial small value (which controls 
the duration of the spinodal stage). In fact the Higgs dynamics has to be 
studied carefully, taking into account, for the Higgs modulus |C|, the expo- 
nential amplification of large wavelength modes due to spinodal instability, 
and the emergence of a homogeneous background out of small quantum per- 
turbations. We will see that these features have a characteristic signature 
on the primordial spectrum : a non-Gaussian spike. 

Let us now start working. The details of our model and its potential are 
given in [23]. However, a numerical simulation of the transition would be 
impossible without a continuous expression for the one-loop corrections around 
the critical value Be, at which symmetry breaking starts. This expression has 
been calculated for toy models of F-term and D-term inflation in [37], and it 
is straightforward to generahze it to our case. The complete effective potential 
is given in the Appendix, together with some remarks on its motivations from 
string theory. We choose a particular arbitrary set of parameters, also given 
in the Appendix. 

In section 2., we recall for non-experts the results usually expected in multiple 
slow-roll inflation. 

In section 3., we treat the problem in the simplest possible way, i.e., neglect- 
ing the Higgs quantum perturbations. Even at this level, we find that the 
amplitude of the large-scale plateau exceeds the analytic prediction. 

In section 4., we discuss the Higgs dynamics at the very beginning of the 
transition, around the unstable equilibrium point. 

In section 5., we go to the next level of complexity for the primordial spectrum 
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calculation, and take into account the Higgs longitudinal perturbations (still 
neglecting the vacuum degeneracy, which is expected here to be irrelevant). 
We use the usual framework of the semi-classical approximation, and carefully 
justify this approach. 



2 Well-known results in multiple slow-roll inflation 



Let us recall some well-know results about quantum perturbations in multiple 
slow-roll inflation. The matter of this section can be found in other papers, 
but we included it for clarifying the results and notations of the next sections. 
We follow essentially Polarski & Starobinsky [38] (see also [39]), and we devote 
a particular attention to a subtle problem (the definition of the V^'s), which 
is of particular importance in our case. 



2.1 Semi- classical equations 



In absence of a phase transition, there are some well-known techniques for the 
calculation of the power spectrum. Usually, the slow-rolling fields (pi^ i=i,..,N, 
are decomposed into a classical zero-mode plus quantum perturbations : 



f (i k 



The field perturbation S(pi, and the metric perturbation (in the longitudinal 
gauge) $ can be expanded over an A'"-dimensional basis of annihilation opera- 
tors aj(k), j=i N, satisfying canonical commutation relations [aj(k), aj/(k)] = 

ScpiiK X, t) = ^ e*^(5(^,,-(^, t)dj{k) + e-'^^-S^^k, t)a]{k). (2) 



This expansion defines a set of mode functions Sipij and $j (from now on, we 
omit all arguments k and t). The basis is chosen in such way that at some 
initial time, the perturbation 5ipi lies along Sj. Commutation relations for the 
fields yield the Wronskian normalization conditions for the mode functions, 
valid at any time : 
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E (^A- = E ^? (3) 

j=i,..,N ^"'^ "^P i=i,..,Ar 

(we are using the reduced Planck mass : mp = (SttG)"^/^ = 2.4 x lO^^GeV). 
The background equations of motion are : 



dV 

with H^ = :^(y + l E*.') W 

(of course, even if not written exphcitly, the potential and its derivatives are 
evaluated at {(pi, ..,(pn))- They imply the useful relation H — — 2^^«'^i- 
The perturbation equations of motion read : 



d'^v\ ■ dV d^V 

^j + H^j = 7^ E (5) 

^"^P i=l,..,N 

Initial conditions for the field mode functions inside the Hubble length (when 
k 3> aH) are given by WKB solutions of equations (5) (with $j neglected), 
normalized with (3) : 

S^,,^S,,^e-^I'^''\ (6) 



Equivalent conditions for the metric mode functions are given at the next 
order in aH/k : 

= i ^ e-'Iadt (7) 



The primordial spectrum is defined as the power spectrum of evaluated 
outside the Hubble radius, during matter domination : 

k^ <0| |$|' |0>=A;3 I^T- (8) 



In fact, in order to obtain the shape of the primordial spectrum, it is not nec- 
essary to follow the evolution of $ until matter domination, but until the end 
of inflation. Indeed, at this time, all the modes observable today are outside 
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the Hubble radius, and are dominated by the so-called "growing adiabatic so- 
lution", which evolution is independent of k. In other words, the primordial 
spectrum J2j changes only by an overall normalization factor between 
the end of inflation and matter domination. 



2.2 Analytic solutions 



We summarize here the analytic solutions for multiple slow-roll inflation with 
uncoupled inflaton fields. It is useful to note that in any case (with or without 
slow-roll), and for each value of j, two solutions of eqs.(5) have an analytic 
expression outside the Hubble radius (when k <^ aH) : 



5^pij _ 1 
(pi a 



t 



Cj J a dt' — Cjj , 



H } . ,\ ~ H 



^.■ = C,|l--yarft' (9) 







The mode with coefficient Cj is the growing adiabatic mode, the one with Cj 
the decaying adiabatic mode. Other solutions (2A^ — 2 isocurvaturc modes for 
each j) have no generic expressions. There is also a freedom to add to $j an 
arbitrary decaying mode V/a. If one looks only for slowly- varying solutions, 
the system (5) simplifies in : 



(10) 

and admits only N slowly- varying independent solutions (for each j). Now, 
during slow-roll stages, we have : 

In this case, the N slowly-varying solutions can be found analytically [38] : 
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The solution with coefficient Cj is the growing adiabatic mode (as can be 
seen by comparing with eq.(9)) and the modes with coefficient are the 
isocurvature modes. For fixed j, only — 1 out of the coefficients dij are 
independent, and the definition is complete only once an arbitrary constraint 
has been chosen (for instance, dij = 0). The system (12) can be inverted in : 



= 2" ^^ij^ 



E,:^ (13) 

A crucial point, to which we will come back later, is that the quantities Vi 
are not uniquely defined. The condition under which the former solution is 
valid during slow-roll reads : Vi, dVi/dt = (dV/dipi) (f^. So, the Vi can be 
found by integrating dV/difi over time or directly over (^j, but there is still a 
freedom to add a constant term to each Vi. This freedom is equivalent to an 
ambiguity in the splitting between adiabatic and isocurvature modes. Indeed, 
we see from eqs.(12) that adding constants to the ^'s is equivalent to changing 
the definition of Cj and dij. So, to remove this freedom, we must give a more 
precise definition of isocurvature modes during infiation. In fact, isocurvature 
modes can be clearly identified at the end of infiation, from the fact that their 
contribution to $j vanishes [38,40]. Using this as a definition, it is possible to 
calculate the V^'s in a unique manner, provided that the whole evolution of 
background fields is known, till the end of infiation (in contrast with single- 
field infiation, for which Cj depends only on background quantities evaluated 
at t = tfc, the time of Hubble radius crossing). 

The most widely-studied multiple infiationary models (for instance, those with 
a polynomial potential and no coupling terms) share the following property. 
When one field exits slow-roll, its time- derivative, which is decreasing, becomes 
small with respect to the other ones ((^^ — > 0, — > est). So, the simplified 
background equations (11) and the solutions (12) hold continuously at any 
time. Simply, at each transition between two stages, the perturbation equation 
for the field exiting slow- roll becomes irrelevant (the perturbation vanishes), 
but the other solutions remain valid. In this case, we have a simple prescription 
for defining the Vi's.[^ In the last stage of single-field inflation, driven by 
one inflaton only the adiabatic mode contributes to $j, by definition. So, 
using eq. (12), we see that we must have Vi = V , and Vi^i = 0. This fixes 
the integration constants : at any time, we can obtain the Vi by integrating 



^ I am very grateful to D. Polarski and A. Starobinsky for providing this 
prescription. 
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backward on the field trajectory, expressed as (<^i, .., (/^at) = Ti{(p,-^, from a 
point at the end of the last slow- roll stage Tj((^j end); back to an arbitrary 
point : 

Vi = 5u ^(r,((^,end)) + / ^iU^i)) d^i- (14) 

end 



So, at each transition, one potential Vi goes to zero, and in the last stage only 
the growing adiabatic mode remains. The interesting quantity at the end, 
which is = (since at the end of inflation the primordial spectrum 

deflned in eq.(8) is equal to C{kY [H /H'^Y), can be immediately found by 
applying eq.(13) at tk. Indeed, at this time, 8ipij ~ '^uy^' 



On the other hand, for models with violent ("waterfall") transitions, i.e., for 
which non-slow-rolling fields acquire large time-derivatives with respect to V, 
the solution of eqs.(12) appUes separately to different slow- rolling stages. Be- 
tween these stages, the background and perturbation evolution can be quite 
complicated. This makes it very difficult to define properly the V^'s, and to dis- 
tinguish between the adiabatic mode and the isocurvature modes (excepted, 
of course, during the last slow-rolling stage). So, the standard result of eq.(15) 
may be difficult to employ, due to the impossibility of defining the V^'s ana- 
lytically. The model considered in this paper is an illustration of this problem, 
which, to our knowledge, had not been addressed previously. 



3 Spectrum neglecting Higgs perturbations 



We will now try to calculate the primordial power spectrum numerically. For 
a pioneering work on this kind of numerical simulation, see ref. [7], in which 
several cases of multiple field infiation are treated. 

The simplest way to study the phase transition is to neglect completely the 
Higgs perturbations, and deal only with the inflaton perturbations {SAj, 5Bj)j=i^2- 
In other words, we consider C as an exactly homogeneous field : C = Ce''^, 
and we solve the following equations : 



A + 3HA+^^0, idem for S,C, 
oA 
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V + 



•-2 --2 

A +B 



■-2 
C 



5 A. + ?,H5A. + + 



dV 



^ dA ^ OAdB ^' 



dV 



OAdB 



(16) 



2m|, 



(AM,- + 555,). 



Note that the small coupling term d^V/dAdB arises only at the one-loop order, 
and when \B\ < Be- At the beginning of the phase transition, C is sitting at 
the origin, which becomes an unstable equilibrium point at the critical e-fold 
number = N^. So, we have to "push" it away, and introduce an initial value 
at Nc, invoking quantum fluctuations as a physical justification. In contrast 
with the results of section 5, the simplified calculation of this section does not 
depend much on this initial value, provided that it is small enough and does 
not introduce a discontinuity. We take however for C the expectation value of 
the coarse-grained quantum fiuctuations, that will be calculated in the next 
section : 



c =< \c, 



quantum \ 



> 



1/2, 



ON i 



(17) 



After Ac, (5, C) rapidly evolve towards the equilibrium point (0, Ceq). In fact, 
it is useful to distinguish five stages : 

(1) first slow-roll stage : A and B are in slow- roll, (7 = 0. 

(2) spinodal stage : while C grows from the initial value given above to a 
critical value Cspinodai{B), the effective mass square is negative : this 
is spinodal instability. This stage will have a characteristic signature only 
in section 5, when C perturbations will be taken into account. In our 
model this stage lasts ~ 0.15 e-folds. 

(3) quasi-static transition : B and C roll down to the minimum. This 
evolution is fast with respect to the Universe expansion (~ 0.8 e-folds 
in our model) and {B, C) are far from being in slow-roll. However, we 
call this transition quasi-static, because for the parameters chosen here 
(in particular (3 ~ 10~^), C remains close to the valley of its local min- 
ima,0 and performs small oscillations around the elliptical trajectory 

= Cl^{\ - B^/Bl) on which dV/dC = 0. 

^ so, we do not deal with the case in which C undergoes a chaotic stage with large 
oscillations. This would arise with 10-2 < /? < 1. In this case, we don't have a 
robust justification for the expression of one-loop corrections [37], and everything 
becomes more comphcated. 
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(4) second stage with oscillations : only A is in slow-roll, but the fields 
{B, C) are not stabilized in their minimum and their oscillations affect the 

background evolution. In our model we find that \C\ <^ \A\ at any time, 
but, if we average over one oscillation, \B\ ^ \A\ during approximately 
15 e-folds. This stage is inflationary | <^ H"^), but differs from a usual 
slow-roll stage since the evolution of H is driven by the oscillating field 
B, not the inflaton field A. 

(5) second stage with slow-roll : at some point, \B\ <ti \A\ and we go 
back to ordinary single-field slow-roll infiation, during approximately 35 
e-folds. 

The slow-roll analytic prediction for the primordial spectrum is given by eq. 
(15), and at first (tree-level) order, the definition of {Va, Vb) resulting from 
eq.(14) and from the potential (eq.(46)) is : 



2(^1 + ^1) 



(18) 



V, 



B ■ 



B^ 



-l3'BC'dB with C^ = max{Ci(l-^),0} 
2 B^ 



during the first stage. 



j3'^B'^Clq{l — ^) during the transition, 

during the second stage. 



(19) 



Note that at any time V — Va-\-Vb- This is due to the fact that the definition 
(14) implies V = Va + Vb + Vc, and if C remains in the valley of minima, 
Vc^J^dCc^O. 



Let us now compare the slow-roll predictions with the results of the numerical 
simulation. We will first discuss in details the results for two modes : on with 
k < a{tc)H{tc) and one with k > a{tc)H(tc). 



3.1 Large wavelength results 



We integrate the equations for a mode that crosses the Hubble radius around 
N — —53, while — —50 (the origin for N has been chosen so that the last 
infiationary stage ends at A'" = 0). The results for |$i| and |$2| are shown 
on the first plot of fig.l. Since the final power spectrum will depend on the 
final value of which is much bigger than |$2| during the last stage, it is 
important to understand the evolution of this quantity, and to compare with 
slow-roll predictions. 
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In first approximation, the evolution of $i can be reproduced analytically. 
Indeed, we check that the perturbation 5-Bi, though amplified exponentially 
during the transition, remains a negligible source term for $i and SAi, except 
during a short period (two e-folds) during which it is comparable with other 
source terms. So, at first order, we can just neglect 6Bi, and deduce the 
evolution of $i and 6Ai from a system of two differential equations instead of 
three.Q. Moreover, we look only for the slowly-varying solution, which obeys 
to : 



There is an exact solution (normalized to \SAi\ = -^^^ at t = t^), that repro- 
duces fairly well the numerical solution : 



oAi = -. ^= exp / T-rzdt, 

-2 t -2 

A 1 r A 

$1 = -. ^ exp / —:r--dt. (21) 

•-2 --2 --2 

We must keep in mind that between and some time ti, one has i? +C > A 
(this is the 15 e-fold "second stage with oscillations"). On the other hand, 

•-2 

during the "second stage with slow-roll", we have 2mpH = —A , and eqs.(21) 
read : 



Hit,) V{h) f B + C' 
(^1 = o ATT^flv/, N exp / T-n-dt. (22) 



Now, since the contribution to the last integral is made essentially for tc < t < 

-2 -2 

ti with H ^ ~ ^2m^ ' ^^^^ term with the exponential can be approximated 
by H^{ti)/H^{tk), and then : 

- ^fa) (23) 



note that this simplification is appropriate for finding and 5Ai, but the 
exact adiabatic solution (9) does not satisfy tfie reduced system (two equations). 
This shows that the solution for $i and 5Ai must be a mixture of adiabatic and 
isocurvature perturbations. 
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This approached result matches exactly the slow-roll prediction (15), because 
V{ti) is nothing else but Va_. In fact, we find that the numerical result (which 

•-2 ^2 ^2 

takes into account the exact relation —2mpH = A + B + C , and the small 
correction due to the fact that 5Bi cannot be neglected during two e-folds) 
exceeds the slow-roll prediction by a factor 4.3. 

To conclude on the evolution of $i, note that it is possible in principle to isolate 
numerically the adiabatic and the isocurvature contribution. We subtract to 
$1 a solution of the type of eq.(9), and also, a solution V/a. We tune (Ci, Ci, V) 
in order to remove any behavior proportional to 1/a or H/a, and to cancel 
$1 during the second slow- roll stage. The result is unique and stands for 
the isocurvature mode. The second plot of Fig. 1 shows the result of this 
decomposition. At the end of this formal exercise, $i appears as the result 
of an almost exact cancellation during the first stage between adiabatic and 
isocurvature contributions, which are of opposite sign. We see that both modes 
are excited during the transition, but the isocurvature mode decays like 
during the second stage and vanishes. 

On the other hand, the solution for j = 2 differs completely from the standard 
slow-roll prediction. We expect from eqs.(15, 18, 19) that C2 ~ 3Ci. It is in fact 
10^ times smaller. Again, this result shows that the usual definition of (Va, ^s) 
is not valid during and before the slow-roll interruption. This appears clearly 
when we separate numerically the adiabatic and isocurvature contribution 
(third plot of fig.l). We see that during the first stage, the isocurvature mode 
completely dominates. 

In conclusion, even if C = (C^ -t-Cl)^/^ differs from the slow- roll prediction only 
by a factor of order one, we have proved that the usual results of cqs.(14, 15) 
cannot be employed in our model. The numerical results can be matched with 
eqs.(13, 15) by tuning Va and Vb during the first stage, but we don't know 
how to predict these values analytically : this calculation seems to be difficult, 
and a numerical simulation is probably unavoidable. 



3.2 Small wavelength result 

We now consider the modes that cross the Hubble radius well after the tran- 
sition. For the smallest wavelengths, this happens during the last single-field 
slow-roll stage, and the usual results of eqs.(15, 18, 19) are automatically vahd. 
In our model, these modes correspond to very small wavelengths today, not 
observable on cosmic scales. 

We are more interested in modes crossing the Hubble radius during the oscil- 
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latory periodQ when the evolution of H is driven by i? ^ A . The evolution 
of ($i,$2) for such a mode, with Hubble crossing at ^ —44 (six e-folds 
after N^)., is shown on the last plot of fig.l. Since we are now in single-field 
inflation, $i and $2 should be pure adiabatic modes; we check that this is 
the case, by fitting ($1, $2) with the expression of eq.(9) (for N > —44). The 
coefficient Ci matches exactly the slow-roll analytic prediction, and C2 ^ Ci, 
as expected for a non-slow-rolling field perturbation. 

This can be explained easily. During the transition, $1 remains a slowly- 
varying solution, because the effective mass of SAi is dominated by k'^/o?'. 
So, at Hubble crossing, $1 is in the slowly- varying adiabatic mode, and Ci 
is exactly the same as in the slow- roll prediction (i.e. eq.(15) with i = 1). 
On the other hand, |$2| is strongly affected by the transition, because the 
effective mass of 5B2 becomes suddenly much larger than fc^/a^. At Hubble 
crossing, |$2| is in the decaying adiabatic mode, proportional to a^^, and sta- 
bilize with a value much smaller than |$i|, as expected for a non slow-rolling 
field perturbation in single-field infiation. 



3. 3 Primordial power spectrum 



We plot on fig.2 the power spectrum of adiabatic fiuctuations, k^^'^C{k), or 
equivalently (up to a constant), A;^/^$(/c). It is a step-like spectrum with su- 
perimposed oscillations on small scales, quite similar to the analytic spectrum 
of Starobinsky [18], but more smooth (the amplitude of the first oscillations 
with respect to the amplitude of the step is smaller). On intermediate and 
small scales, it is also close to the spectrum of [9] (from double "polynomial" 
infiation), with the notable difference that on large scale we have an approxi- 
mately flat plateau, instead of a logarithmic dependence on k. 

For our choice of primordial parameters, the step has got an amplitude p ~ 10, 
and spans over one decade in k space. As we saw before, this spectrum arises 
solely from Ci, i.e., from metric perturbations calculated with initial condition 
of the type {5 A ^ 0, 55 = 0). 



Since the oscillations of B produce oscillations in the effective mass of C, one may 
expect that SC modes could undergo a stage of parametric resonance, that would 
show up in resonant metric amplification [43]. This is not the case because the field 
C is very heavy, so that the relative amplitude of the mass oscillations is kept very 
small. 
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4 Dynamics of the Higgs field 



Until now, we have treated the Higgs field C as a simple classical homogeneous 
background field, and put by hand an initial condition ai N — Nf., in order to 
push the field away from unstable equilibrium. We will now enter into more 
details concerning the Higgs quantum fluctuations, in order to justify and 
calculate the initial condition at Nc (by averaging over quantum fluctuations), 
and to prepare the work of the following section, in which we will consider the 
effect of the Higgs quantum perturbations on the primordial power spectrum. 
We rewrite the potential (Appendix, eq.(46)) under the more suggestive form : 

V = V{A, B) + \(5\B' - Bl)C' + (24) 
When B < Be, the potential is a usual "Mexican hat" with respect to C. 



4-1 Fluctuations during the first stage of inflation 



We decompose the field C as in eqs.(l,2), using now a three-dimensional op- 
erator basis dj^ i=i,2,3- Since C is complex, we would need in principle a four- 
dimensional basis, in order to quantify separately the real and imaginary part 
Ci and C2. However, in this section, the U{1) symmetry will still be preserved, 
and the mode functions are identical for Ci and C2. So, for concision, we in- 
troduce only one degree of freedom j = 3, describing simultaneously both 
directions in the complex plane. 

The Wronskian condition (3) applies also to the mode function SCj. During 
the first stage of inflation, C = 0. So, in addition to eqs.(16) (now taken for 
J = 1, 2, 3), we have the following perturbation equation : 

SCj + 3HSC, + [^-^ + 0^1 = 0, (25) 



with a WKB solution valid both within and outside the Hubble radius : 



with 



'2m 

4 ■ 



(26) 



During the first stage of infiation, at the first perturbative order considered 
here, $3 vanishes, and so do the non-diagonal mode functions (^^3, 5Bz, 
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SCj^z)- is decoupled from other perturbations and its evolution depends 
mainly on the mass d'^V/dC'^ = — B'^) that goes to zero at = Nc- 



4-2 Formation of inhomogeneities just after Nc 



When > N^, the mass of C becomes negative, causing an exponential 
amplification of the mode functions, known as spinodal instability. Let us 
restrict the analysis to the very short stage after during which SC^ remain 
small enough for the linear equation (25) to be still valid.[^ As we can check 
a posteriori, this stage lasts ~ 0.15 e-folds. If we re-scale SC3 to x = a^^'^SC^, 
we see that eq. (25) can be rewritten : 



We will solve this equation under a few approximations, and compare our 
result with an exact numerical simulation. Since we are interested in a short 
period after Nc, with typically N — Nc < 0.2, it is appropriate to linearize the 
A^-dependence of the effective mass-term in equation (27) : 



-fi\N-Nc), /i' 



d{d^V/dC^} 



dN 



> 0, 



(28) 



'1-2{N-Nc)) 



(in the particular model we are studying, n can be easily found from the slow- 



/3^B dV 
3^2 dB 



/32 2\n2{g\U+9UB) 



roll condition for 5 : /i^ - - ^ 

short stage, H is approximately constant and N — Nc 
approximations, equation (27) reads : 



Also, during this 
H{t~tc). Under these 



X - fjl H {t - tk) X = 0, 



(29) 



with : 




^il=f? + 2'-, tk = tc + -i^A'---.H''\- (30) 



^ this is the case if (0| |/ (2^^)3/2 i^^fc Pi 0)^^^ ^ Ceg : then, the cubic term in 
dV/dC is negligible with respect to the linear term, and f^(/ j^^^Ck) ~ 
0(0) X J^6C,. 
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The modes x start growing exponentially at t = tk (for k = 0, t = tc— j^, but 
smaller wavelengths start growing later). The correctly normalized solution of 
equations (29) and (25) is given, up to an arbitrary complex phase, by : 



X-- 



(31) 



6Cs = a-'/\, 



where H^^fz] is a Hankel function of the first kind. Indeed, for t < tk, this 
solution has an asymptotic expression which is identical to equation (26). For 
t ~ tk, one has : 



6C. 



-3/2 _ 



1/3 



3r(|) \fikHy\ 



(32) 



and for t > tfc in very good approximation : 



, expf|^fc/7V2(t_t )3/2 

^2/Xfci71/2(t_t,)l/2 



Is it possible to define an effective background from the coarse-graining of long 
wavelength modes? The field C, coarse-grained over a patch of size / ~ 2Tia/k, 
is a stochastic classicalf^ quantity. Its real and imaginary part obey to a 
Gaussian distribution (at least if all modes are in the vacuum state), with 
variance <C'^>^/'^ computed from : 



The effective background can be considered as classical when large wavelengths 
modes have very large expectation values. Then, observers measuring the value and 
momentum of C cannot feel the non-commutative operator structure of the field, 
and the quantum field behaves as a classical stochastic field. More precisely, it is 
know [41] that the modes can be approximated by classical stochastic quantities 
when their expectation values are much bigger than the minimal expectation value 
set by Heisenberg uncertainty principle : 



(^1 y(k)pt(k) + C.C. 1^-) > (*| [y(k),pt(k)] 1^- 



1. 



(34) 



This condition on operators can be translated in terms of mode functions and read 
for the vacuum state : 



\5C6C* + 5C*5C\ > -i{5C5C* - 5C*5C). (35) 

For eq.(26) (before Nc) this means H ^ mc, which does not hold for the non slow- 
rolling field C\ but when the exponential amplification starts, the condition above 
is rapidly satisfied for spinodal modes. 
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<C'>^ II ^^(0|5C(k,x,t)5Ct(k',x,t)|0) 







27r 



k^\6Csfdk, (36) 







where keg is the coarse-graining cut-off. When t < tc, we know that scales 
outside and around the Hubble radius have a ^-independent amplitude (see 
eq.(26)). So, <C^> scales like k^^ (intuitively, all modes have the same am- 
plitude, but small wavelength have a bigger statistical weight, proportional to 
k"^). On the other hand, when t > tc, large wavelengths start being amplified 
earlier; so, in the coarse-graining integral, there is a competition between the 
term |(5C3p and the term fc^, and a scale will be privileged. More precisely, 
let us use the asymptotic expression (33). The function k'^lSC^l'^ peaks around 
the scalef^ 

= «^ ( V^^) ' ^^^^ 



and only modes in the range O.lkmax < k < "ik^ax, which are within the 
Hubble radius, contribute to the coarse-graining integral (36). A numerical 
simulation confirms this result for 0.1 < N — < 0.2 (see fig.3). This means 
that just after Nc, the field C has an inhomogeneous structure, and can be 
seen as effectively homogeneous only in regions of typical size A <C a/kmax- 

Inside each patch, we can define an effective homogeneous background from the 
coarse-graining of large wavelength modes, and use the standard semi-classical 
approximation. The constant phase in a given patch can be chosen randomly, 
while the squared modulus \C\ obeys a probability distribution that can be 
calculated exactly. Indeed, with any coarse-graining cut-off keg ~ 'ikmax, we 
can compute the root mean square of the real and imaginary part of C : 

Since Ci and C2 are Gaussian stochastic numbers with variance < >^/^, 
then |Cp obeys a distribution and \C\ has a mean value ^ y^3/2 <C^ >^/'^ 

and a variance ^ y^l/2 < >^/^ (the numerical factors can be found in 
tables; see, for instance, eq. (26.4.34) in [42]). We see explicitly, by comparing 
with eqs.(32, 33), that < >^/^ is a solution of the zero-mode equation 

to find this it is appropriate to use the approximation fe^lCsp cx 
Pexp(4/3/zi/V2(i_i^__g^)3/2)_ 
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with initial condition : < >V2 (tj = {6TiT{l))-\3fi'^Hy/^. Actually, the 
numerical simulation shows that the exact result exceeds this estimate by a 
factor 2.3, so that the correct initial condition is : 

<C^>^/^ {Q = 0.12 {^^Hf/\ (39) 

So, we really have in each patch an effective homogeneous background, obeying 
for > 0.1 to the equation C + + dV/dC = 0, and with an initial 
distribution of probability. 

On the other hand, on scales k < kmax, the background is essentially inho- 
mogeneous, and due to strong mode-mode coupling, the linear semi-classical 
approximation breaks soon after Nc- 



4-3 Later evolution and the homogeneous phase approximation 

So, it is not possible to find an exact solution for the evolution of C dur- 
ing spinodal instability and at later time. We should recall here that out-of- 
equilibrium phase transitions and spinodal instability have been widely studied 
in condensed matter physics, and also in the context of "spinodal inflation", 
using either the large-N limit [31,32] (which is well-defined, but appropriate 
for the symmetry breaking of U{N), N —>■ oo), or just a Hartree-Fock ap- 
proximation [33] (which is more difficult to interpret, and cannot account so 
far for the metric back- reaction). In such studies, the typical scale of homo- 
geneity for the coarse-grained background is larger than the Hubble radius 
{kmax < o,{tc)H{tc)), while the observable primordial spectrum arises from 
modes with k > a{tc)H{tc)- So, it is possible to separate completely the modes 
contributing to the "effective zero-mode assembly" , and those contributing to 
the primordial spectrum. In contrast, in our case, we need to follow modes 
on scales k with k < a{tc)H{tc) < kmax- So, at first sight, it seems that our 
problem is much more complicated than spinodal inflation, due to the inho- 
mogeneity of the background. In fact, it is not, because we have an infiaton 
field A driving inflation even during the transition, and we will see that a 
perturbative semi-classical approach can still be employed.P^ 

The following observations will suggest an approximation under which we can 
continue the primordial spectrum calculation : 



Another difference with respect to spinodal inflation is that we will never coarse- 
grain or average modes far outside the Hubble radius. This is just because the 
spinodal field is not slow-rolling, and then, modes are exponentially amplified also 
on sub-Hubble scales. 
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• since the causal horizon during the first inflationary stage is much bigger 
than the Hubble radius, {A, B) are homogeneous on the largest scales ob- 
servable today, Be is reached everywhere at the same time, and the phase 
transition is triggered coherently. So, the modulus \C\ should be approxi- 
mately homogeneous, even on the largest scales : we should be able to cast 
it into a zero-mode plus small perturbations. 

• the phase inhomogeneities should not be much relevant for the primordial 
spectrum calculation. The existence of a negative mass, leading to exponen- 
tial amplification, only arise in the longitudinal direction \C\. The degree 
of freedom associated with the complex phase is a Goldstone boson; it is 
eaten up by the gauge field, which becomes massive. Starting from zero, 
this mass soon exceeds the Hubble parameter : so, there cannot be large 
quantum fiuctuations arising from this sector. The existence of phase in- 
homogeneities is also linked to the formation of cosmic strings at the very 
beginning of spontaneous symmetry breaking, separated at least by a char- 
acteristic length A > a/kmax- These strings could eventually contribute to 
the formation of cosmological perturbations today, but as we said in the 
introduction, we do not include such mechanism in this study. 

So, for our purpose, it seems reasonable to neglect the phase inhomogeneities, 
forget the vacuum degeneracy, and assume a homogeneous phase 6'(x, t) = 9. 
In the following section, we will recompute the primordial power spectrum, 
taking now into account the Higgs longitudinal quantum perturbations, and 
we will find that these perturbations have a characteristic signature on the 
primordial power spectrum. 



5 Spectrum with Higgs longitudinal perturbations 

We go on neglecting the vacuum degeneracy. We saw in section 4.2 that just 
after the transition, the modulus \C\ has got a mean value ^^3/2 < |Cp >"'^/^ 

and a variance y^l/2 < |Cp >^/^ that we calculated explicitly. Assuming a 
homogeneous phase is almost equivalent to identifying this modulus with a 
real scalar field, canonically quantized, and decomposed, following the usual 
semi-classical approximation, into a background field C = y^3/2 < |Cp >^/^, 
plus quantum fiuctuations. The quantum mode functions should be matched 
with the quantities SCj that we already studied around Nc, divided by -\/2. 
With these prescriptions, the mean value and the variance of the modulus 
previously studied in the last section coincide exactly with those of the real 
scalar field introduced in this paragraph.P^ 

at least, when {N — Nc) > 0.1 (i.e., when < |Cp >i/2 obeys the zero-mode 
equation, see eq.(38)), and (A^ — Nc) < 0.2 (i.e., when the equation for SCj remains 
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This matching is the key point of our study. It is not exact, because, in a clas- 
sical stochastic sense, it is based only on the first two momenta of the distri- 
bution : we artificially "gaussianize" the fluctuations. However, the matching 
seems appropriate for estimating the amplitude of the quantum fluctuations 
of the modulus |C|, even if any information on a possible non-gaussianity is 
lost. Since in this approach we introduce a (physically justifled) zero-mode, 
we expect that the growth of spinodal modes will slow-down and terminate 
rapidly, and that the perturbative semi-classical approach will remain valid. 
This will be justifled explicitly in section 5.3. 

In summary, we simulate the usual background equations for [A, B, C, H) 
and the following perturbation equations : 



5Aj + ?,H5Aj + 
5Bj + ZHSBj + 



5Cj + ?,H5Cj + 




dV d'^V 



dB ^ 



KA . 

dAdB ^ dBdC 



dV 

4(7$ - 2 $ - 

^ ' dC ' dBdC 



5Bj, 



2ml 



{ASAj + BSBj + CSCj 



When we start following a mode 6Cj well inside the Hubble radius, we employ 
the initial condition of eq.(26) divided by a factor a/2. At N = Nc, we put by 
hand a non zero- value for C, given by eq.(39) multiplied by a factor \/3/2. 



5.1 Large wavelength results 



The evolution of (|$i|, |$2|, l^sl) for a mode crossing the Hubble radius at 
~ —53 (i.e. during the oscillatory stage) is shown on the upper plot of 
fig.4. By subtracting the previous solution (fig.l) to the new solution (fig.4), 
we check that the evolution of |$i| is unchanged, except during a few e-folds 
after the beginning of the transition. The difference is due to the small (one- 
loop) coupling term qj^^ : when 5Bi is amplified, it excites |$i|, but later 
on, this mode decays and doesn't leave an observable signature. On the other 
hand, the negative effective mass of 5Cj and the large (tree-level) coupling 
between 5Bj and 5Cj generate a considerable amplification of |$2| and |$3|. 
At the end, when the isocurvature modes and the decaying adiabatic mode 
can be neglected, the remaining growing adiabatic mode is smaller than the 
one for |$i|, but only by one order of magnitude. 

linear, see eq.(25)). 
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So, on large scale, what is the observable difference between these results, and 
those of section 3 ? First, since initial conditions for 6B2 are proportional to 
{a\/k)~^, we expect the large-scale power spectrum for A;^/^|$2| to be, as usual, 
approximately scale-invariant. In other words, the contribution found for one 
mode will be the same for other large-scale modes. So, the large-scale power 
spectrum will still be a flat plateau, but possibly with a small deviation from 
gaussianity, caused by the small contribution of |<l'2|- As we said before, our 
method does not allow a quantitative evaluation of this deviation, and we leave 
this for further studies. Second, since initial conditions for SC3 are independent 
of the scale, and since the factor /c^/a^ is always sub-dominant in the effective 
mass of SCs on large wavelengths, |$3| should be scale- independent on these 
wavelengths. So, the contribution of /c^/^|$3| to the power spectrum k^^"^ < 
|$p >^/2 T^ill be proportional to /c^^^, and will peak for scales crossing the 
Hubble length during the transition. We therefore expect an observable non- 
Gaussian spike, whose shape and amplitude will be calculated in section 5.4. 



5.2 Small wavelength result 

The evolution of (|$i|, |$2|, l^sl), for a mode crossing the Hubble radius at 
~ —43 (i.e. during the oscillatory stage), is shown on the lower plot of fig.4. 
and I $2 1 ai'G identical to previous results of section 3.2. |$3| has the same 
behavior as |<l>2|, but is even smaller. So, for small wavelengths, including the 
Higgs modulus perturbation makes no difference at all. 



5.3 Consistency check of the semi-classical approximation 

Before giving the primordial spectrum of perturbations, we will check here 
the consistency of our approach, and the validity of the semi-classical approx- 
imation in our case. Strictly speaking, the semi-classical approximation holds 
whenever for each field, the Klein-Gordon equation 

f)V 

g-'/%{g'/'g'^''dM^, t)) + ^ (<^(x, t)) = (41) 



can be casted into a background equation for the zero-mode, plus independent 
linear equations for field and metric mode functions. The first term, with space- 
time derivatives, is linear provided that metric perturbations remain small : 

^^^^ 
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The second term, with the potential derivative, is hnear in two cases 



(1) obviously, when the perturbations are small with respect to the zero- 
mode : 

<0| I / (^f)3/2 (5(/?|^|0>^/^< ip. Then, we can write : 

- (^(x,t)) c -(^) + ^(^) X y ^5^. (43) 

(2) also, when both the zero-mode and the perturbation expectation value 
are small enough for non-quadratic terms in the potential to be neglected. 
Then, we have the following linearization : 

dV d'^V 

5^(^(x,^))= 5^X¥^(x,t), (44) 

and d^V/dip"^ is a constant. 

We have seen that around = A^c + 0.15, in the homogeneous phase approx- 
imation, we can define a zero-mode and a set of perturbations. Since all these 
quantities are very small, we are in the second case, and the semi-classical 
approximation is valid. But what happens later? We know from figure 4 that 
one e-fold after Nc, the perturbations of C, B and $ reach a maximum; if, 
at that time, the expectation values of the perturbations exceed the homoge- 
neous background fields (or, in the case of metric perturbation, exceed one), 
then our approach is simply inconsistent. 

To address this question, we first study the evolution of the power spectra 
of SCj, 6Bj and for A^c < -/V" < A^c + 1- The results for 5Cj are shown 
on three-dimensional plots, in fig.5. This simulation is just the continuation 
of the one performed in section 4.2, and at low A^ we recognize the curves 
of fig. 3. On the upper plot, we show J2j < I^C'jP >i logarithmic scale, 
vs. k and A^. We see that large scales (spinodal modes) are exponentially 
amplified during the spinodal stage, and then perform coherent oscillations, 
driven by the background evolution. On the lower plot, we show k^J2j < 
|5Cjp >, in linear scale, vs. k and A^. We observe that the main contribution 
to the averaging integral (eq.(36)) always arise from the same region in k 
space, centered around the approximately constant wavenumber kmax- This 
observation is technically important, because it justifies the use of finite limits 
of integration in the averaging integral. In particular, we see that we don't 
have to worry with ultra-violet divergences. Provided that the upper limit of 
integration kcut-off is chosen in a finite range (namely, around the largest k 
values shown of fig.5), then < |Cp > does not depend on this cut-off. Similar 
results are found for 6Bj and . 
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Then, by integrating on k, we find the expectation values < |Cp >^/^, < 
|^|2 >i/2^ ^ |<|,|2 >i/2^ g^]^^ p^Q^ them with respect to C, B and 1 (see fig.6). 

It appears that for N c:^ Nc + 0.20, when C approaches the valley of min- 
ima (on the figure, we recognize this valley from the small oscillations around 
it), i.e., when dV/dC becomes non-linear, the quantum perturbations are sub- 
dominant by one order of magnitude. When C reaches Ceq, the quantum expec- 
tation value is smaller by three orders of magnitude. After, we know from fig. 4 
that the field will become more and more homogeneous. For 5, we see that the 
perturbations are always sub-dominant (when B performs its first oscillation, 
it goes through zero, but the perturbations are still quite small and we remain 
in the linear regime). For the maximal value is < |$p >i/2r^ 10"^. So, the 
semi-classical approach holds at any time in the range < N < 1, and 
therefore, until the end of infiation. 

In conclusion, in this model, the semi-classical approximation can be employed 
during the whole evolution, because there exists a time (A^ ~ A^^^ -|- 0.15) at 
which the perturbations have grown sufficiently to create an effective zero- 
mode (obeying the usual zero-mode equation), but not too much, so that we 
are still in the linear region of the potential derivative. It is at this timeP^ that 
we have done the matching between the complex field (with no zero-mode) and 
a real scalar field (with a zero-mode) standing for its modulus. If such a time 
did not exist, complicated non- linear approaches and further approximations 
would have to be employed. This would be the case if the second infiaton 
field A was not supplying a constant potential energy, and supporting the 
Universe expansion. Indeed, in that case, the evolution of the scale-factor and 
of the fields would become completely stochastic, and it would be a hopeless 
challenge to keep track of modes crossing the Hubble length before and during 
the transition. 

5.4 Primordial spectrum including Higgs longitudinal perturbations 

We show our final result on figure 7, for the parameters used previously (mid- 
dle). On the other plots, we show some primordial spectra corresponding to 
different values of (3. The contribution of |$i| is unchanged with respect to fig. 
2, but we now have an additional contribution from $2 and $3 on intermediate 
and large scales. Let us summarize the characteristics of the power spectrum : 

• On small scales (large k), we have an approximately flat plateau, exactly 
identical to the analytic predictions. Its amplitude can be found in our previous 
paper ([23], eqs. (6,7)), as a function of the Fayet-Illiopoulos terms {^a, C,b), 

this time can be clearly identified on the upper plot of fig.6 : without the factor 
v/372 (resp. 7172) in C (resp. < |C7|2 >V2)^ the curves would be tangent around 






iV = + 
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of the coupling constants (qa, Qb) and of the duration of the second stage AiV 
(different values of AA^ correspond to different choices for the inflaton fields 
initial values). The plateau has got a tilt = 1 — 1/ AN ~ 0.98, like in usual 
(single-field) inflationary models, when the potential derivative is provided by 
logarithmic one-loop corrections. 

• On intermediate scales (the value of these scales is given by AN), we find 
a spike, as expected from sections 5.1 and 5.2 (since A;'^/^|$3| is proportional 
to A;^/^ for small k, and strongly suppressed for large k). This spike may be 
non-Gaussian, for reasons that we already explained. Clearly, it is beyond 
the scope of this paper to quantify the non-gaussianity, because the semi- 
classical approximation is not appropriate for extracting the statistics of the 
modes. The amplitude and width of this spike with respect to large scales 
depend essentially on /3. Indeed, P is proportional to the Higgs mass, and 
then, controls the duration of the transition : for low /3 values, the transition 
is long, exponential amplification is enhanced, and the spike is higher. Its width 
is roughly of one decade in k space (but, since A;^/^|$3| is proportional to k^^"^ 
on large scales, small spikes are automatically more narrow). The shape of the 
spike around its maximum is slightly non-trivial, reflecting the first oscillations 
of the field B around zero during the transition. It is reasonable to state that if 
the primordial power spectrum really possesses a spike (in the range in which 
this is still an open possibility z.e., 0.1 x lO^^/iMpc"^ < A: < 5 x lO^^/iMpc^^), 
then its amplitude is already constrained by observations to be smaller than 
~10 (here, we define the amplitude as the ratio between the maximum value of 
A;^/^|$| and the one, for instance, of the large-scale plateau). This corresponds 
to the condition /3 > 10~^. 

• On large scales (small k), we also have an approximately flat platcaii, hut as 
we saw in section 3.1, we don't know how to predict analytically the amplitude 
and tilt of this plateau, due to the violent slow-roll interruption during the 
transition. For our choice of parameters {^a, Cb, 9a, 9b, ^N), there is a ratio of 
4.3 between the amplitude obtained numerically for and the analytic 
prediction (eq.(13) of [23], or eq.(23) of this paper). This ratio is found to 
be approximately independent of /?. The large scale tilt is very close to one. 
Finally, metric perturbations on these scales are composed of a (Gaussian) 
component and of a (possibly non-Gaussian) component $2, generated, like 
$3, by the exponential amplification of perturbations during the transition. So, 
the relative contribution of $2 depends, like the spike, on /?. For f3 > 10~^, we 
find that $2 is sub-dominant, as can be seen in the figures, but not negligible 
(for P — 10~^, |<I>i|/|$2| = 3.4). We also investigate the dependence of the 
large-scale plateau on ^a, Cb- Since small-scale amplitude depends on ^a — Cb 
(see eq.(7) of [23]), we vary ^a + ^b while keeping ^a — ^b fixed. We find that 
the amplitude on large and intermediate scales increases with C,a + Cb- On 
fig.8, ^A + Cb has been divided by 2.8, and other parameters are the same as 
in fig.7 (upper plot). In this case, the ratio between the two plateaus is 1.5. 
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6 Conclusion 



We calculated numerically the primordial power spectrum for a particular 
model of double supersymmetric inflation, called double D-term inflation [23]. 
Following the usual semi-classical approach, we integrated for each indepen- 
dent mode function the equations of motion of field and scalar metric pertur- 
bations. Some technical difficulties (attached to any model of double hybrid 
inflation) arise between the two inflationary stages. Indeed, the first stage ends 
with a rapid second-order phase transition, taking place during inflation, and 
corresponding to the spontaneous symmetry breaking of a U{1) gauge sym- 
metry, triggered by a Higgs field C. This transition is followed by another 
inflationary stage, equivalent to usual hybrid inflation. We were especially 
interested in the behavior of modes exiting the Hubble radius during the tran- 
sition, since they are expected to produce non-trivial features in the primordial 
power spectrum. 

After spontaneous symmetry breaking, it is impossible to follow exactly the 
mode functions of the complex field C. Indeed, on the largest observable wave- 
lengths, the complex phase of C is completely inhomogeneous, at least at the 
beginning of the phase transition; so, on these scales, the usual perturbative 
semi-classical approach is broken by mode-mode coupling. However, we ar- 
gued that only longitudinal Higgs perturbations, i.e., perturbations of \C\, 
can contribute to the primordial power spectrum. Indeed, these perturbations 
undergo a stage of spinodal instability, during which large wavelength modes 
are exponentially amplified, while the perturbations of the complex phase and 
of the gauge field (which combine into a massive gauge field, following the 
usual Higgs mechanism), are those of a non-slow-rolling field with large posi- 
tive mass square. 

So, we performed a calculation in which longitudinal Higgs perturbations are 
taken into account, with the modulus \C\ treated as an ordinary real scalar 
field, described by the semi-classical equations. We showed that at the begin- 
ning of the transition, a zero- mode for |C| emerges from quantum fiuctuations 
of C's real and imaginary parts. Since the final result depends crucially on this 
initial zero-mode, we calculated it precisely in section 4.2. Then, by compar- 
ing the background fields and the averaged quantum perturbations, we showed 
that the semi-classical approach holds at any time. However, during one e-fold 
or so after the beginning of the transition, it is not far from its limits of validity 
(with typically a factor 10^ between background quantities and perturbations). 
This fact is a first hint that the contribution to the power spectrum resulting 
from longitudinal Higgs perturbations could be significantly non-Gaussian. 
The second hints comes from the fact that at the beginning of the transition, 
we identified a modulus with a real scalar field. In this process, we matched 
carefully the mean value and expectation value of both quantities, but infor- 
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mation on higher momenta was lost. It is beyond the scope of this paper to 
evaluate the non-gaussianity. 

Our final power spectrum, shown on fig. 7 for different parameter values, can 
be divided in three regions : 

(1) a small-scale Gaussian plateau, whose amplitude and tilt {ns — 0.98) can 
be calculated analytically. 

(2) a possibly non-Gaussian spike, with an amplitude depending on the su- 
perpotential coupling parameter /3. For (3 ~ 10~^, the spike induces vari- 
ations at most of order one in the primordial spectrum. 

(3) a large-scale plateau, for which we don't have an exact analytic prediction. 
For /3 ~ 10"'^, deviations from gaussianity, if any, are expected to be 
small. 

For a precise comparison of our model with CMB and large-scale-structure 
(LSS) data, one needs in principle, not only the primordial power spectrum 
of scalar adiabatic perturbations, but also the spectrum of primordial gravi- 
tational waves, and the contribution of cosmic strings formed during the two 
phase transitions (the one considered here, and the one occurring at the end 
of inflation). However, as we said in [23], the tensor-to-scalar ratio T/S (for 
the first temperature anisotropy multipoles) depends essentially on the gauge 
coupling constants (as in single-field D-term infiation). Both qa and cjb should 
be at most of order 10^^ (if not, infiation takes place with infiaton values of 
the order of the Planck mass, and supersymmetry/supergravity theories are 
not expected to be valid; sec the Appendix for the consistency of this require- 
ment with string theories). With this order of magnitude for the coupling 
constants, one can show easily that T/S is negligible in double D-tcrm in- 
fiation. On the other hand, cosmic strings can contribute significantly to the 
matter power spectrum and to CMB anisotropics, especially for D-term infia- 
tion, as shown in [27] (sec also [28]). Since, in current studies, perturbations 
from strings and from the primordial spectrum are just added in quadra- 
ture, it is perfectly reasonable, at least in a first step, to calculate only the 
primordial spectrum. The ratio of strings-to-infiation large-scale temperature 
anisotropics, Rsi — /C^, is difficult to estimate; theoretical predictions 
give a ratio around 3 [27] or 4 [28] , but the authors of this last reference argue 
that there are many uncertainties in these calculations, so that smaller values 
can still be considered seriously. With a mixture of infiationary and string 
perturbations, our model could still be interesting : the features would appear 
nevertheless, but smoothed by the cosmic string contribution. 

In this work, we did not compare our results with CMB and LSS observations. 
Indeed, current data cannot distinguish with precision the kind of small fea- 
tures that we obtain, at least if they are located on scales 0.1 x 10"^/iMpc~^ < 
A; < 5 X 10~^/iMpc~^. However, as we said in the introduction, comparisons 
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with observations of typical BSI spectra (with a step or a spike) have already 
been performed. The conclusions are currently quite encouraging (see for in- 
stance [17,15]) and provide strong motivations for studying BSI inflationary 
models, until a final statement is made by future redshift surveys, such as the 
2-degree Field projectf^ or the Sloan Digital Sky Surveyp^ combined with 
forthcoming balloon and satellite CMB experiments. 

An important qualitative aspect of our results is the possible deviation from 
gaussianity on intermediate and large scales. This may be related with the 
evidence for a non-gaussianity in the COBE data [28], but further work is 
needed on this point, since, on the one hand, we don't have a precise theoret- 
ical prediction of non-gaussianity in this model, and on the other hand, the 
conclusion of [28] involves large uncertainties. 

Finally, we should stress that experimental evidence for a BSI primordial 
power spectrum would be a very positive breakthrough for cosmology. On 
the one hand, the introduction of one or two additional inflationary parame- 
ters, associated with the shape of the spike or the step, would not compromise 
cosmological parameter extraction, as indicated by [44], because the effect of 
these parameters would be orthogonal to the other's. On the other hand, exper- 
imental data would encode more information, and provide additional exciting 
constraints on inflationary models, fundamental parameters, and high-energy 
physics. 
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Appendix 



The model. Double D-term inflation [23] is based on the superpotential 
W = aAA_^.A_ + PBB^B^, and on the definition of two gauge groups {U{1)a, 
U{1)b), with corresponding gauge coupling constants {qa, Qb) and Fayet- 
Iliopoulos terms (^a, ^b)- We choose the charges to be (0,0) for A, B, (±1, 0) 
for A-t, and (±1,±1) for B± (other choices can also lead to successful models). 
During inflation, we need only to follow the fields A, B and B_. Moreover, 
the complex phases of A and B remain constant during inflation (except for 
very peculiar initial conditions that we don't consider here). So, we deal with 
four canonically normalized real fields (^4, B, Ci, C2) defined as : 

A = V2\Al B = V2\B\, Ci + iC2 = C =V2B_. (45) 



The scalar potential reads : 



+ 



(46) 



The last term is the one-loop correction. Following [37], (a^, m^) are taken 
from the following table : 



i 


at 




1 


2 


la'A' + g\{U-m') 


2 


2 




3 


2 




4 


2 


\P\2B^ + {Cf) 


5(» 


2 


¥'B^ - glU - gUn 


(<) 


1 


2{glU + g%^B-'^P'B') 


6(» 


2 





«) 


3 




7 


-4 


\a^A^ 


8 


-4 


{\(3^B^ + \{I3^ + 2g\ + 2gl){Cf + Va) 


9 


-4 


(1/3^5^ + |(/3^ + 2g\ + 2gl){CY - 



30 



Here, (rcsp. means "when \B\ > B" (resp. "when \B\ < B^^), and 



Be 



^ '^(9aU+9b ^^ ^j^g table, we have used the notations : 



(47) 



\ i9l + 9l) 

A = 1 {P^B^ + iP^ + 2g\ + 2gl){Cyy - \(5\g\ + gl){C)\ (48) 

The one-loop corrections are automatically continuous in = Bf.. The renor- 
malization scale A must be chosen around A ~ {3 Be and, as in [37], we optimize 
this choice numerically by requiring the continuity of the potential derivative. 

Choice of parameters. The dependence of the primordial spectrum on the 
parameters is discussed in section 5.4. Before this section, we choose some 
arbitrary parameters (Ca, 5'_b) inside the allowed region defined in [23], 

for which the power spectrum has got the same order of magnitude as the one 
indicated by COBE, and the amplitude of the step (between the large-scale 
and the small-scale plateau) is of order one : 



/a = 3 X 10-^mp, JiB = 4.2 X IQ-^mp, gA = gB = lO"'. (49) 



The value of a is completely irrelevant for our study (it is relevant only at 
the very end of inflation). On the other hand, the choice of (3 is crucial. Until 
section 5.3, we take /3 = 10~^. As mentioned in section 3, footnote 6, values 
of (3 greater than 10^^ would render the problem very complicated; we would 
not have a robust justification for the expression of one- loop corrections, and 
the transition would consist in chaotic oscillations in (S, C) space. 

String motivations. Double D-term inflation is a consistent model from the 

point of view of supersymmetry, when the Fayet-Iliopoulos terms are put by 
hand from the beginning. Here, we briefly address the issue of consistency 
with string theory. It turns out that in heterotic string theory, one can always 
redefine the gauge groups in order to end up with only one anomalous [/(I), 
and a single associated Fayet-Iliopoulos term. So, in this framework, it is im- 
possible to obtain two ^ terms at low energy. However, it is well-known that 
even the simplest models of single D-term infiation are hardly compatible with 
heterotic string theory [45,46], mainly because the Fayet-Iliopoulos term gen- 
erated by the Green- Schwartz mechanism is typically too large to account for 
COBE data. It has been noted by Halyo [47] that the situation is much better 
in type I string orbifolds, or type IIB orientifolds. Then, the Fayet-Iliopoulos 
terms are not generated at the one-loop order [48], but at tree level, and de- 
pend on the vacuum expectation value of some moduli [49] corresponding to 
the "blowing-up modes" of the underlying orbifold; since there is some free- 
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dom in fixing the moduli, the value of ^ required by COBE does not appear 
as unnatural in this framework (also, in these theories, there is no unification 
of all coupling constants as in heterotic string theory, and one can envisage 
scenarios where the coupling constant g of the anomalous gauge symmetries is 
slightly lowered with respect to the Standard Model couphng constant [47,50], 
so that the infiaton takes values significantly smaller than the Planck mass 
during inflation, and supergravity, or even global supersymmetry can be em- 
ployed). But, nicely, another unusual feature in type I string orbifolds and 
orientifolds of type IIB is the coexistence of several anomalous [/(l)'s, with 
associated Fayet-Ihopoulos terms [49]. So, our model may naturally arise from 
these theories. We are very grateful to Stephane Lavignac for having kindly 
provided information on these points. 
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Fig. 1. Evolution of the metric perturbations, when Higgs perturbations are ne- 
glected : (top) |$i| and |<I>2| for a long wavelength mode; (middle left) splitting of 
|$i| in adiabatic and isocurvature contributions for this mode; (middle right) same 
for |<&2|; (bottom) |$i| and \^2\ for a small wavelength mode. The transition starts 
at Nc — —50, while the last inflationary stage ends at ~ 0. The long wavelength 
mode exits the Hubble radius at = —53, the small wavelength mode at = —44. 
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Fig. 2. The primordial power spectrum of adiabatic fluctuations, k^^'^C{k), or equiv- 
alently, /c^/^|$(A;)|, when the Higgs perturbations are neglected (with logarithmic 
scale and arbitrary units). It is a step-like spectrum with superimposed oscillations 
on intermediate scales. 
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Fig. 3. The left column shows the Higgs perturbation power spectrum : log |(5C3p 
versus log/c, for different times chosen around the beginning of the transition, i.e, 
when the effective mass of C becomes negative. From top to bottom, we have 
N — Nc = —0.025,0.075,0.175. The vertical scale is the same in the three plots. 
We show on the horizontal axis the particular values corresponding to A; = aH and 
k = a|m|, where m? = d'^V/dC'^. We see that spinodal modes, with k < a\m\, are 
amplified. On the right column, we plot the integrand of the coarse-graining inte- 
gral : /c^|(5C3p. The vertical scale increases from top to bottom. It appears clearly 
that after a fraction of e-fold, only modes with H < k/a < \m\ contribute to 
this integral, which indicates the emergence of an inhomogeneous background, with 
first-order homogeneity recovered on smaller scales (higher k). On the last plot, we 
also show the scale kmax calculated analytically from eq.(37). 
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Fig. 4. We now take into account the Higgs longitudinal perturbations, and show the 
evolution of the metric perturbations (|^i|, |^2|) l^aD) ^or (top) a long wavelength 
mode, crossing the Hubble radius at = —53, and (bottom) a small wavelength 
mode, crossing at = —43. 
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Fig. 5. Before calculating the expectation value of quantum perturbations, and 
compare it with background quantities, it is useful to follow the various power 
spectra, plotted here for 6C. The upper plot, |5Cp (logarithmic scale) versus k 
and A^, shows that spinodal modes are exponentially amplified, and then oscillate 
coherently. The second plot, A;^|5Cp (linear scale) versus k and A'^, shows that the 
main contribution to the averaging integral arises always from the same modes, with 
k ~ kmax] this plot shows how to define the limits of integration in the averaging 
integral. 
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Fig. 6. Consistency check of the semi-classical approximation. The expectation value 
of the quantum perturbations < \(p\^ >i/2 ig generally much smaller than the back- 
ground fields (than one in the case of metric perturbations). When it is not the 
case, both quantities are very small and the linear approximation is still valid. The 
three plots correspond to the case of C, B, and $. The case of A is not shown 
because we know from the beginning that this field remains extremely homoge- 
neous. The evolution is shown only until N = Nc + 1, which is the time of maximal 
inhomogeneity, as can be seen on fig. 4. 
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Fig. 7. The primordial power spectrum of adiabatic fluctuations, k^^'^C{k), or equiv- 
alently, k^^'^^{k), when the Higgs longitudinal perturbations are taken into account 
(in logarithmic scale and arbitrary units). Together with the total power spectrum, 
we plot the three contributions from <^j=i,2,3- For scales crossing the Hubble length 
during the phase transition, <I>3 produces a spike (possibly non-Gaussian). On the 
first three plots, the parameters {Ca,S,b, dA, 9b) are fixed to the values of the ap- 
pendix, and /? = 2 X IQ-^ (top), lO'^ (middle), 0.5 x IQ-^ (bottom). 
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Fig. 8. The primordial power spectrum of adiabatic fluctuations, for the same pa- 
rameters as on fig. 7 (upper plot), excepted lower values of S,a, ^b- On the right, the 
small scale plateau is unchanged, because we kept — £,b fixed. On the left, we get 
less power on large and intermediate scales. 
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